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Background & Motivation 


• When chemical kinetic timescales are approximately equal to flow 
timescales, the chemical composition of a flowfield must be 
determined as part of a simulation procedure. Such flows are in 
chemical nonequilibrium. 
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When chemical kinetic timescales are approximately equal to flow 
timescales, the chemical composition of a flowfield must be 
determined as part of a simulation procedure. Such flows are in 
chemical nonequilibrium. 


The physical models and governing equations for flows in 
thermochemical nonequilibrium have been simulated previously with 
finite difference and finite volume techniques. 

In this work we review the physical models and implement the first 
known SUPG finite element scheme for hypersonic flows in 
thermochemical nonequilibrium. 


T 



• Molecules and atoms can store energy in 
various modes. 


• At hypersonic conditions these modes may not 
be in equilibrium, resulting in thermal 
nonequilibrium. 



March 23, 2011 4 / 25 




Physical Modeling 


Governing Equations 


Governing Equations 

• Extension from a single-species calorically perfect gas to a reacting 
mixture of thermally perfect gases requires species conservation 
equations and additional energy transport mechanisms. 


dp 

dt 

dpu 

~dT 

dpE 

~df 


V ■ {p u) =0 

V ■ ( puu ) = —VP + V • r 

V ■ ( pHu ) = —V • q + V • (tu) 


March 23, 2011 5/25 


Kirk, Bova, Bond 


SUPG FEM for Thermochem. Noneq. Flows 




Physical Modeling 


Governing Equations 


Governing Equations 

• Extension from a single-species calorically perfect gas to a reacting 
mixture of thermally perfect gases requires species conservation 
equations and additional energy transport mechanisms. 

(p s u) = V ■ (pV s Vc s ) + lj s 
( puu ) = —VP + V ■ r 

( ns 

pJ2h s V s Vc s 

S=1 


dps 

dt 

dpu 

~dT 

dpE 

~dt~ 


V ■ 


V ■ 


V ■ 


March 23, 2011 5/25 


Kirk, Bova, Bond 


SUPG FEM for Thermochem. Noneq. Flows 




Physical Modeling 


Governing Equations 


Governing Equations 

• Extension from a single-species calorically perfect gas to a reacting 
mixture of thermally perfect gases requires species conservation 
equations and additional energy transport mechanisms. 

^ + V ■ (p s u) = V • ( pV s Vc s ) + 
at 

dpu 


dt 

dpE 

~dt~ 


+ V • ( puu ) = —VP + V • r 

/ ns 

+ V • ( pHu ) = -V • q + V • (tu) + V ■ I h s V s Vc 


S=1 


• Problem class may also require a multitemperature thermal 
nonequilibrium option. 


dpey 

dt 
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Thermochemistry 


Thermodynamics & Transport Properties 

• Thermochemistry models must be extended for a mixture of 
vibrationally and electronically excited thermally perfect gases. 


Here we have assumed that T trans = T rot = T and T vib = T elec = T v 
• Additional transport property models are required. In this work we 
use 

► species viscosity given by Blottner curve fits, 

► species conductivities determined from an Eucken relation, 

► mixture transport properties computed via Wilke's mixing rule, and 

► mass diffusion currently treated by assuming constant Lewis number. 


e int =e trans + e rot + £ vib + e elec + ^0 


,trans 


ns 


= 5>eT s (T) + E c s ef(T) + 


.trans 


s—1 s—mol 

ns ns 


ns 


E c «ef (T v ) + J2 c s e T C (T v ) + J2 c sh° s 
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Physical Modeling 


Thermochemistry 


Chemical Kinetics 

• We consider r general reactions of the form 


N 2 + X ^ 2N + X 

N 2 + O — NO + N 


• The reactions are of the form 


TZ r 



kf r n 

S = 1 



Ot S r 


where a sr and /3 sr are the stoichiometric coefficients for reactants and products 
• The source terms are then 


nr 

UJs — M s ^ ^ (c^sr fisr') (J^br 
r= 1 
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Physical Modeling 


Thermochemistry 


Kinetic Rates 

• The forward rate coefficients are defined with a modified Arrhenius 
law as a function of some temperature T 

k fr (?) = C fr T' n ' r exp ( -E ar /RT ) 


where the rate constants are determined empirically. 

• The corresponding backward rate coefficient can be found using the 
principle of detailed balance and the equilibrium constant K eq 


K — til 
£q hr 


• In thermal equilibrium T = T. We are currently using Cantera in this 
regime. 

• In thermal nonequilibrium f = T (T, T v ) and typical hackery ensues. 
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Physical Modeling 


Thermochemistry 


Energy Exchange 


Cjy — Q v + ^transfer 
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Physical Modeling 


Thermochemistry 

Energy Exchange 

UJy = Qv “L Qtransfer 

We adopt the Landau-Teller vibrational energy exchange model 


e vib = ps~- 


r vib 


( 1 ) 


where e" b is the species equilibrium vibrational energy and the vibrational 
relaxation time r^ lb is given by Millikan and White 


_vib _ 

' e 


E ns 7i T ( ns 

r= lXr _ M „ 

YTrUXr/T^ XV ° r M r ’ 


Cs 


and 


exp [An- (t~ 1/3 - 0.015^/ 4 ) - 18.42 
= 1.16 x 10 - 3 / 4 / 2 0 4 / 3 , Hsr = 


vib ±_ 

1 sr p 


M s + M r 
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Physical Modeling 


Thermochemistry 


Vibrational Energy Production and Energy Exchange 

Wy = Q v + Qtransfer 




Physical Modeling 


Thermochemistry 


Vibrational Energy Production and Energy Exchange 

Wy = Q v + Qtransfer 

When molecular species are created in the gas at rate 6j s , they contribute 
vibrational/electronic energy at the rate 

Qvs = Us (ef + ef c ) 

so the net vibrational energy production rate is 

ns 

Q„ = ^(ef + ef c ) (2) 

5=1 
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Physical Modeling 


Thermochemistry 


Vibrational Energy Production and Energy Exchange 

Wy = Q v + Qtransfer 

When molecular species are created in the gas at rate 6j s , they contribute 
vibrational/electronic energy at the rate 

Qvs = Us (ef + ef c ) 

so the net vibrational energy production rate is 

ns 

Q„ = ^(ef + ef c ) (2) 

5=1 

Combining (1) and (2) yields the desired net vibrational energy source 
term 

ns ns 

^ = + (ef + ef c ) 

S=1 5=1 
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Finite Element Formulation 


dU dFi dGi ■ 
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Finite Element Formulation 
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Finite Element Formulation 


dU . AC .pv dU d / dU\ ■ 

+ ( Aj + Aj ) ( Kij ' ) + <S 


dxj dxj 


dx n 


Find U satisfying the essential boundary and initial conditions such that 
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for all W in an appropriate function space. 
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Finite Element Formulation 



March 23, 201 1 


12/25 


Kirk, Bova, Bond 


SUPG FEM for Thermochem. Noneq. Flows 




Results 


Inviscid Thermal Nonequilibrium Chemically Reacting Flow 


Inviscid Cylinder 

• Inviscid flow in thermochemical nonequilibrium 

• 5 species air (N 2 , 0 2 , NO, N, O) 

• 5 reaction model with Park 1990 rates 

cN 2) oo = 0.78,cO 2iOO = 0.22 
Uoo = 5, 500 m/sec 
Poo = 10“ 2 - 10 -4 kg/m 3 
Too = 250 K = Tv, oc 

• Landau-Teller vibrational energy relaxation model, with Millikan and 
White vibrational relaxation time 
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Inviscid Thermal Nonequilibrium Chemically Reacting Flow 
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Results 


Inviscid Thermal Nonequilibrium Chemically Reacting Flow 
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Results 


Inviscid Thermal Nonequilibrium Chemically Reacting Flow 



t£ = ^ exp [A sr ( t ~ 1/3 - 0.015 M ^ 4 ) - 18.42] 
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Results 


Viscous Thermal Equilibrium Chemical Reacting Flow 

2D Extended Cylinder 

• Laminar flow in thermal equilibrium 

• No-slip, adiabatic, noncatalytic wall 

• Chemical nonequilibrium, 5 species air (N 2 , 0 2 , NO, N, O) 

• 5 reaction model with Park 1990 rates 

cN 2j00 = 0.78 ,cO 2i oo = 0.22 
t/oc = 0- 731 m /sec 
Poo = 6.81 X 10“ 4 kg/m 3 
= 265 K 

• Blottner/Wilke/Eucken with constant Lewis number Le = 1.4 for 
transport properties 

• Mesh, iterative convergence 

• FIN-S/DPLR comparison 

• Weak & Strong Scaling 

x 
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Results 


Viscous Thermal Equilibrium Chemical Reacting Flow 
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Results 


Viscous Thermal Equilibrium Chemical Reacting Flow 


Mesh Convergence 
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Relative Transient Residual, \\U/\t\ 


Results 


Viscous Thermal Equilibrium Chemical Reacting Flow 


Iterative Convergence 



Time Step 
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Time Step Size, At, (seconds) 




Results 


Viscous Thermal Equilibrium Chemical Reacting Flow 


Speedup 





y (rn) 


Results 


Viscous Thermal Equilibrium Chemical Reacting Flow 


Code-to-Code Comparison - 

Stagnation Line 



■ i I i 7 ■ i I i i i i l 

0.00 0.05 0.10 0.15 

X (m) 



x(m) 


March 23, 2011 22/25 


Kirk, Bova, Bond 


SUPG FEM for Thermochem. Noneq. Flows 


o z z 




y (m) 


Results 


Viscous Thermal Equilibrium Chemical Reacting Flow 
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Near-term Effort 


Additional Focus Areas 
O Physics Modeling 

► Weakly Ionized Flows 

► Surface Catalycity 

► Additional Boundary Conditions 

© Ablation coupling 
© Adjoints 

► Sensitivity analysis 

► Adaptivity 



Primal & Dual Error 



Combined Qol Error 
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Near-term Effort 


Thank you! 
Questions? 
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